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ABSTRACT 

The specific incidence of radio flares appears to be significantly larger than that of 
the prompt optical emission. This abundance, coupled with the reverse shock inter- 
pretation suggests that radio flares add a unique probe on the physics of GRB shocks. 
Motivated thus, we estimate the strength of the reverse shock expected for bursts 
in which multi- wavelength observations have allowed the physical parameters of the 
forward shock to be determined. We use all 6 bursts (980519, 990123, 990510, 991208, 
991216, 000418) which are found to be adiabatic and thus predicted to have a strong 
reverse shock. We aim to constrain the hydrodynamic evolution of the reverse shock 
and the initial bulk Lorentz factor - which we found to be between 10 2 and 10 3 and 
well above the lower limits derived from the requirement that gamma-ray bursts be 
optically thin to high-energy photons. In half of the cases we improve the description 
of the early afterglow lightcurves by adding a contribution from the reverse shock. 
Modelling of this early emission provides the opportunity to investigate the imme- 
diate surroundings of the burst. For 991216 and 991208, the expected 1/r 2 density 
structure for a stellar wind is not compatible with the early afterglow lightcurves. 
Considering the radial range relevant to these GRBs, we discuss the conditions under 
which the inclusion of a wind termination shock may resolve the absence of a 1/r 2 
density profile. 
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1 INTRODUCTION 

The prompt and extremely bright optical flash in GRB 
990123 (Akerlof et al. 1999) was accompanied by a strong ra- 
dio flare (Kulkarni et al. 1999). Peaking at ~1 day, the flare 
was unlike usual radio afterglows which rise to maximum on 
a timescale of weeks or even months. The simplest interpre- 
tation of this radio flare is that it arises from the reverse 
shock component of the external shock (Sari & Piran 1999, 
hereinafter SP99). The reverse shock (hereinafter RS) prop- 
agates into the adiabatically cooled particles of the coasting 
ejecta, thereby shocking the shell material and producing a 
prompt optical flash. The RS emission then weakens rapidly 
and shifts to lower energies, eventually crossing the observed 
radio band. Indeed, it has been shown that the radio coun- 
terpart to GRB 990123 is compatible with estimates derived 
from scaling down the prompt optical emission to the epoch 
of the radio afterglow observations (Kulkarni et al. 1999). 

It turns out that there is a broad range of model param- 
eters for which a strong optical flash is expected to precede 
the main GRB afterglow (Soderberg & Ramirez-Ruiz 2002; 
hereinafter SR02) . Such flashes are, however, difficult to de- 
tect in the current observing modes (Akerlof et al. 2000, 



Paczynski 2001). On the other hand, the specific incidence 
of radio flares of 1:4 (Djorgovski et al. 2001) appears to be 
significantly larger than that of the prompt optical emis- 
sion obtained by ROTSE (Kehoe et al. 2001) or LOTIS 
(Williams et al. 2000). This abundance, together with the re- 
verse shock interpretation, suggests that the radio flare phe- 
nomenon has the potential to shed new light on the physics 
of GRBs. Motivated by this interpretation, we have esti- 
mated the strength of the radio flare expected from bursts in 
which broadband observations have been able to constrain 
the physical parameters of the forward shock (hereinafter 
FS) emission (Panaitescu & Kumar 2001b; Panaitescu & 
Kumar 2002, hereinafter PK02). We use all 6 bursts (980519, 
990123, 990510, 991208, 991216, 000418) which are found to 
be adiabatic (i.e. v m < v c ) for reasonable assumptions about 
To and thus predicted to have a strong RS. Radiative bursts 
(e.g. 970508) are not selected since the reverse shock compo- 
nent is significantly quenched in this regime, and therefore 
the observations in these cases are dominated entirely by 
emission from the forward shock. By evolving prompt optical 
flux estimates to the epoch of radio afterglow observations, 
we are able to discern whether a contribution from the RS 
could have been detected. We find that for ^ of the bursts 
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for which the cooling frequency is larger than the typical 
synchrotron frequency, the predicted RS emission (which is 
generally found to lie in the mildly relativistic temperature 
regime) combined with that of the FS produces an improved 
fit to the radio afterglow data. As we shall discuss, the strong 
dependence of the peak time of this radio flare on the bulk 
Lorentz factor F provides a way to measure this elusive pa- 
rameter. Accepting the conclusion that the radio flare arises 
from the RS, we then place constraints on the velocity of 
source expansion and hence on the density and profile of 
the medium in the immediate vicinity of six bursts (980519, 
990123, 990510, 991208, 991216, 000418). Fhese constraints 
provide environmental diagnostics for GRBs and are com- 
pared with the main types of environments considered for 
afterglows. We assume Ho = 65 kms -1 Mpc -1 , Qm = 0.3, 
and Qa = 0.7. 



2 THE REVERSE SHOCK: OPTICAL FLASH 
TO RADIO FLARE 

Observations of the optical flash associated with GRB 
990123 confirmed earlier predictions of prompt emission 
from a RS (Meszaros & Rees 1999, hereinafter MR99; SP99). 
Although it contains a comparable amount of internal en- 
ergy to that of the FS, the RS bears a temperature which is 
significantly lower (typically by a factor of F) . As it crosses 
the shell, the RS produces a single short burst of emission 
which is predicted to peak in the optical band with reason- 
able assumptions for the burst physical parameters (MR99). 
The time of peak emission, t pC ak, strongly depends on the 
initial bulk Lorentz factor of the burst. For high To the peak 
time is comparable to the duration of the burst while for low 
To it typically occurs at later times. More specifically, the 
peak time is defined as 



peak 



max[A/c, tdec], 



(1) 



where A is the width of the shell. Here, tdec, is the time at 
which the inertia of the swept-up matter significantly decel- 
erates the shell ejecta 



tde 



3E 
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1/3 



(2) 



where E is the isotropic energy of the burst and no char- 
acterises the density of the external medium (SP99). A can 
be inferred directly from the observed burst duration (tdur) 
by noting that A = ctdur/(l + z) and assuming the shell 
does not undergo significant spreading (Piran 1999). At a 
given time (e.g. r pea k), the broadband RS spectrum can be 
described by the ordering of the three synchrotron break fre- 
quencies: the self-absorption frequency v a , the cooling fre- 
quency v c and the characteristic synchrotron frequency v m . 
For spectra with v a < v m < v c , the RS is adiabatic while 
for v a < v c < v m the RS is radiative, and so the electrons 
are cooled quickly. Although a radiative RS could demon- 
strate a peak flux of comparable brightness to an adiabatic 
blast wave, the emission would be rapidly quenched making 
it exceedingly difficult to detect. 



2.1 Relativistic versus subrelativistic 

The RS spectral break frequencies are easily calculated by 
comparing them to those of the FS (MR99; SP99; Panaitescu 
& Kumar 2000) . By assuming equality of velocity and pres- 
sure across the contact discontinuity separating the shocks, 
it is possible to define the properties of the reverse shocked 
region in terms of no and To (Blandford & McKee 1976; 
hereinafter BM). However, unlike the FS, the RS is not al- 
ways relativistic (Sari & Piran 1995). Shells satisfying 



E 



1/6 



x A 



-l/2 F -4/3 



>> 1 



(3) 



are thin and are likely to have a Newtonian RS, which is 
typically too weak to decelerate the shell effectively. On the 
contrary, for thick shells (£ << 1), the RS is relativistic and 
thus successfully described by BM solution. For an adiabatic 
blast wave peaking at v m , the thin shell spectral energy 
equations are given by 

r.9 2 



v m = 5.8 x 10 a el_ 1 ea / l 2 no / rg(l + 2 )" 1/4 Hz (4) 



Fl, max = 4.2 x 10 °D 2S -e B _ 2 



1 x 10 8 e 2 _ 1 e^ 2 _ 2 n V2 r 2 (l + ; ? ) 



1/2 „l/2 r ?f1 _1_ ^-V4] 

i 2 e^ /2 _ 2 F;5onJ /2 ro(l + «) 3 / 8 Jy, (5) 
while for the thick shell case we have 

1/2 Hz (6) 

F,, max = 6.0Di a e^ 2 _ a ^ 4 nS /4 ro J C /4 (1 + ( 7 ) 
where F Vilnetx is the spectral peak flux, e e and ts are the 
equipartition functions for the electrons (e) and for the mag- 
netic field (B) respectively, and D is the luminosity distance 
of the burst. Here we adopt the convention Q = 1Q X Q x for 
expressing the physical parameters. It is important to note 
that in the thin shell regime -F„ iinax scales directly with To 
while for the thick shell regime the relation is inverted. 

2.2 Light curves of the reverse shock emission 

Unlike the synchrotron spectrum, the afterglow light curve 
at a fixed frequency strongly depends on the hydrodynam- 
ics of the relativistic shell, which determines the temporal 
evolution of the break frequencies. In a mildly relativistic 
RS (£ > 1), the temperature of the shocked material is 
non-relativistic and so the late-time evolution of the ejecta 
no longer follows the BM solution. Furthermore, the Sedov- 
Taylor solution is not applicable due to the the relativistic 
bulk Lorentz factor of the fluid (Kobayashi & Sari 2000; 
hereinafter KS00). The hydrodynamic evolution of a mildly 
relativistic RS therefore lies in a regime for which there are 
no analytic solutions available. In order to constrain the evo- 
lution of T in this regime it is common to assume (MR99; 
KS00) 



r oc R~ 9 with (3/2 < g < 7/2) 



(8) 



The limits on g reflect two evolutionary pathways of the 
pressure. Adiabatic expansion (p oc p 4 ^ 3 ) is assumed for 
g — 3/2 while g — 7/2 corresponds to the case of pres- 
sure equilibrium, i.e. when the pressure of the FS matches 
that of the ejecta. Using the relation t oc R/Y c, one can 
obtain the scaling of V in terms of the observer-frame time, 
which is given by 

roct- 9/(1+2s) . (9) 
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The pressure and density then scale as 



p oc i 



-8(3+g)/7(l+2g) 



p OC t 



-6(3+g)/7(l+2 9 ) 



(10) 



Application of these scalings to a slow cooling RS gives the 
spectral evolution as a function of g. For an adiabatic blast 
wave, it can be shown that v m oc Tp 5 / 2 p~ 2 and F Vm oc Tp 1 ^ 2 
(KS00). In terms of the observer-frame time, these scalings 
are then 



x t -3(8+5 9 )/7(l+2 9 ) 



F„ m oct-( 12+ll9 ^ 1+2B \ (11) 



Following the peak optical emission, with the passage of 
time, the light from the RS shifts rapidly to lower frequen- 
cies. After the RS peak frequency v m crosses the observing 
frequency it decays rapidly. The evolution of the radio flux 
for frequencies above v m is typically found by assuming a 
distribution of injected electrons with power law of index 
p above a minimum Lorentz factor 7*. The corresponding 
spectral flux at a given frequency above v m is then given by 
F v « Fy m (z//i/ m ) _( - p_1 ^ 2 ; with F v oc v 2 at low frequencies 
{v < Va) and F„ oc v 1 ^' i bridging the two regimes. Together 
with the scalings of equation 11, the above expression gives 
the evolution of the radio flare immediately after the v m 
crossing: 



F x T -(7+24p+15p9)/14(l+2 9 ) y > 



(12) 



It should be noted that for a typical spectral index p = 2.5, 
the flux decay index varies in a relatively narrow range 
(w 0.4) between the theoretical limits of g = [1/2,3/2]. 
Since the flux decay index is a monotonic function of g it is 
not very sensitive to it, and consequently the decay shape 
is more strongly affected by the value of To than of g. As 
demonstrated in Section 3, the quality of GRB radio obser- 
vations is not yet sufficient to place robust constraints on g, 
and thus the hydrodynamic evolution of the RS is currently 
not well determined. 



3 CONSTRAINING THE KINEMATICS AND 
HYDRO-DYNAMICAL EVOLUTION OF THE 
REVERSE SHOCK 

Despite ongoing observational searches, GRB 990123 re- 
mains one of few bursts for which a radio flare has been 
observed. The ability to detect flares depends on both the 
strength and peak time of the optical flash as well as the 
subsequent flux decay index. In turn, these observables de- 
pend on the ambiguous dynamics of the reverse shock. By 
constraining the kinematic and hydrodynamic properties of 
the shock, we may enable observational strategies for radio 
flare detection to be improved. Moreover, as we will demon- 
strate, constraints on the physics of the RS may provide 
us with environmental diagnostics of GRBs. In this section, 
we present a method to determine the properties of the RS 
through modelling of the early radio afterglow observations. 
We constrain both the velocity of source expansion To and 
the evolutionary index g of the RS in bursts for which multi- 
wavelength observations have been able to determine the 
physical parameters of the FS. We use all 6 bursts (980519, 
990123, 990510, 991208, 991216, 000418) which arc found to 
be adiabatic (i.e. v m < v c ) for reasonable assumptions about 
To. Assuming that the burst equipartition values for the FS 
are shared by the RS, we apply the formalism developed in 
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Figure 1. GRB 980519: constraints on To and g assuming a 
thin shell. The best fit value is marked by a solid dot and the 
approximate 1 a uncertainty contour is shaded. A RS contribution 
can be minimised for values of To ~ 490 and g ~ 1.5. Thick shell 
equations provide a lower limit of To > 230 (filled arrow) which 
is consistent with To > 250 derived by PK02. This RS is found to 
be mildly relativistic. Afterglow observations are fit better with 
a FS model (lower panel, solid line) than with a FS + RS model 
(lower panel, dashed line). 



the previous section to determine the spectral and tempo- 
ral evolution of the RS. We estimate the strength of the RS 
contribution in the radio band during the epoch of afterglow 
observations. Through combining our radio flare estimates 
with predictions for the FS emission (obtained by PK02), 
we are able to estimate the combination of To and g which 
provides the best % 2 fit to the early radio data. We consider 
both cases where the temperature of the shocked shell is 
relativistic or not. This is particularly important since for 
reasonable assumptions about the velocity of source expan- 
sion we find in all cases that £ ~ 1. For Xfs+rs/Xfs < 1j the 
fit is improved by including emission from the RS. We will 
show that although current RS theory predicts a detectable 
radio flare component for a wide range of To and g param- 
eter values, it appears to be observationally supported for 
only some bursts. Below we discuss the radio emission and 
RS modelling of the individual GRBs. A summary of results 
are listed in Table 1. 



3.1 GRB 980519 

The broadband modelling of the FS emission by PK02 was 
found to overestimate the radio observations at early times. 
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The addition of a RS contribution to the FS emission serves 
only to worsen the fit. Assuming the thin regime solution, 
we find that a minimum flux contribution is expected from 
the RS for To ~ 490 and g ~ 1.5. Thick shell solutions for 
this burst do not converge since the predicted radio flare 
decreases indefinitely with increasing IV Therefore, we can 
only determine a lower limit of To > 230 for this regime. 
These constraints imply a mildly relativistic RS with £ ~ 1. 
In this case, neither the thick nor thin regime are ruled out 
by the observations. For these values of To and g, the radio 
observations are fit with a FS + RS emission model giving a 
similar \ 2 to that of the FS alone. We find a reduced \ 2 = 12 
and a ratio of of the quality of the fits gives Xfs+rs/Xfs = 
1.4 for the thin shell case (see Figure 1). It should be noted 
that the x 2 ratio can in principle decrease for values of To > 
1000 and g = 1.5, however, with a peak time of only a 
few seconds, this solution is clearly inconsistent with the 
constraints given by equation 1. Therefore, by assuming the 
results for Xfs+rs/Xps ~ 1.4 as our best fit, we find this 
solution to be in agreement with the previous estimate by 
PK02 of To > 250. 

Here we comment on the structure of the \ 2 region for 
this burst as a representative of the larger sample included 
within this study. First, it should be noted that the unusual 
shape of the \ 2 region is the result of an intrinsic coupling 
between the fitted parameters, To and g. Although a more 
proper treatment of the uncertainty region would include a 
full Monte Carlo analysis to determine the error region, this 
level of complexity is not warranted by the quality of avail- 
able data. We therefore apply standard \ 2 fitting techniques 
to build two dimensional uncertainty regions and emphasize 
that these represent only approximate confidence regions. 
Current theory dictates that the uncertainty region for To 
has a lower boundary of To ~ 30 due to the lack of photon- 
photon attenuation (Meszaros, Laguna & Rees 1993) and an 
upper boundary of To ~ 10 3 based on the pulse width evo- 
lution of the prompt emission (Ramirez-Ruiz & Fenimorc 
2000). These theoretical bounds should be adopted as the 
overall uncertainty region for each of the bursts within this 
sample. Secondly, it should be noted that the sharp inflec- 
tion points traced by the confidence region are an artifact 
of the discontinuity imposed by the transition from a non- 
relativitsic to a relativistic regime. The tendency for the 
best-fit value of g to favour the parameter space extrema 
demonstrates that Fo is the dominant parameter in this 
model and thus the Lorentz factor drives the spectral evo- 
lution of the RS emission within this parameter space (as 
discussed in Section 2.2). 

3.2 GRB 990123 

Early radio observations of GRB 990123 revealed a flare 
rising to maximum after ~ 1 day and then rapidly fading 
away (Kulkarni et al. 1999). The strong dimming of the radio 
emission after ~ 2 days excluded a FS origin. The simplest 
interpretation was that it arises from the RS (Kulkarni et al. 
1999; Panaitescu & Kumar 2000). As expected, the FS mod- 
elling by PK02 significantly underestimates the radio emis- 
sion at the time of the flare. By including a RS contribution, 
a FS + RS emission model improves the fit by a factor of 
~ 100. Best fit values are obtained for T ~ 850 (r ~ 240) 
and g ~ 3.4 (g ~ 3.3) for the thin (thick) shell case (Fig- 
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Figure 2. GRB 990123: Constraints on To and g assuming a 
thin shell regime. Best fit values arc marked by a solid dot and 
approximate 1 a uncertainty contours are shaded. Radio flare 
parameters which produce the best fit to the data I\) ~ 850 and 
g ~ 3.4 for the thin shell case. Best fit FS + RS light curve 
models arc plotted (dashed line) within the lower panel in order 
to compare with FS model alone (solid line). The FS+RS model 
provides a significantly better fit to the radio emission. (Note the 
discrepancy associated the first point is due to the effects of self- 
absorption.) The To constraints are consistent with the estimate 
of To > 300 (upper arrow) derived by PK02 and the lower limit 
of To > 180 (lower arrow) by LS01. Note the mildly relativistic 
nature of the thin shell regimes, with £ ~ 1. 



ures 2 and 3) . Both regimes predict a mildly relativistic RS 
with £ values of ~ 1 and ~ 6 for the thin and thick shell 
solutions, respectively. Despite the differences between the 
regimes, both solutions show a minimum at Xfs+rs = 0.3 
which gives a ratio Xfs+rs/Xfs ~ 0.01, thus demonstrating 
a remarkable improvement to the FS fit. Note that the in- 
clusion of self-absorption effects, which are not included in 
the FS model of PK01 or the FS+RS model presented here, 
may enable a better description of the radio observations 
before one day (the reader is referred to SP99 for this mat- 
ter). It should be commented that within the la uncertainty 
region of the thin regime solution, there is a secondary min- 
imum which occurs for lower Lorentz factors, near Fo ~ 400 
and g ~ 3.5. This solution corresponds to a peak time of 
200 seconds which is clearly inconsistent with the observed 
optical light-curve which peaks at ~ 50 seconds (Akerlof et 
al. 1999). Based on the observed behaviour of the optical 
flash, which is successfully described by a mildly relativis- 
tic shell (see Kobayashi 2000), we thus conclude that the 
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Figure 3. GRB 990123: Constraints on To and g assuming a 
thick shell regime. Best fit values are marked by a solid dot and 
approximate 1 <r uncertainty contours are shaded. Radio flare 
parameters which produce the best fit to the data To ~ 240 and 
g ~ 3.3 for the thick shell case. Best fit FS + RS light curve 
models are plotted (dashed line) within the lower panel in order 
to compare with FS model alone (solid line). It should be noted 
that the thin shell equations produce better fits to both radio and 
optical emission (see Figure 2). The To constraints are consistent 
with the estimate of To > 300 (upper arrow) derived by PK02 
and the lower limit of To > 180 (lower arrow) by LS01. Note the 
mildly relativistic nature of the thick shell regimes, with £ ~ 6. 

thin shell solution is the regime which better describes the 
RS behaviour. We find the values of To consistent with pre- 
vious estimates for this burst: To ~ 1200 (Wang, Dai & 
Lu, 2000), T = 1400 ± 700 (Panaitescu & Kumar, 2001a), 
To « 900 ± 100 (SR02), To > 180 (Lithwick & Sari 2001, 
hereinafter LS01), and r > 300 (PK02). Estimates for the 
dimensionless thickness parameter, £, are similarly in agree- 
ment and include: £ ~ 0.7 (KS00) and £ = 0.99 (SR02). 
Furthermore, the hydrodynamical evolution has previously 
been characterised by g = 2.2 (KS00) which is consistent 
within our approximate 1 a uncertainty region. 

3.3 GRB 990510 

As in the case of GRB 980519, the broadband modelling of 
the FS emission of this burst overestimates the radio obser- 
vations at early times (PK02) . We find that a minimum flux 
contribution from the RS can be attained from the thin shell 
equations by setting To ~ 270 and g ~ 1.5. The thick shell 
solution provides only a constraint of To > 230. This RS is 



Figure 4. GRB 990510: Constraints on To and g assuming a 
thin shell. The best fit value is marked by a solid dot and the 
approximate 1 a uncertainty contour is shaded. The radio flare is 
minimised for values of To ~ 270 and g ~ 1.5 with £ ~ 2. Thick 
shell equations provide a lower limit of To > 230 (filled arrow) 
which is consistent with To > 140 derived by PK02 (upper unfilled 
arrow) and Tq > 79 by LS01 (lower unfilled arrow). Afterglow 
observations are better fit with a FS model (lower panel, solid 
line) than with a FS + RS model (lower panel, dashed line). 

found to be mildly relativistic with £ ~ 1. A FS + RS model 
gives a best fit value with x 2 — 8.5 and Xfs+rs/Xfs ~ 8-0 
indicating that the standard FS fit is still preferred (see Fig- 
ure 4). A better \ 2 ratio (~ 1) can be obtained for values 
of To > 2000 and g — 1.5. We disregard this solution, how- 
ever, as it corresponds to i pC ak < 1 s. Note that these results 
are in agreement with the estimate of Po > 140 derived by 
PK02 and the lower limit of r > 79 by LS01. 

3.4 GRB 991208 

Broadband fits to the FS of GRB 991208 included two sets 
of radio lights curves observed at 8.5 GHz and 15 GHz. We 
calculate radio flare predictions for both observing frequen- 
cies. At 15 GHz, the FS prediction underestimates the data 
(PK02). By adding a RS contribution, we find the ratio 
Xfs+rs/Xfs can De improved by a factor of ~ 2 for both 
thin and thick solutions. Best FS + RS fits are produced 
using To ~ 130 and g = 3.5 for the thick shell equations 
which give Xfs+rs ~ 1-5 and a ratio of Xfs+rs/Xfs ~ 0.4. 
The thin shell solution is similar in quality, with best val- 
ues of To ~ 110 g ~ 1.9, giving Xfs+rs ~ 1-9 an d a ratio 
of Xfs+rs /Xfs ~ 0.5. Both regimes predict a mildly rela- 
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Figure 5. GRB 991208: Constraints on To and g assuming a 
thin shell. The best fit value is marked by a solid dot and the 
approximate 1 a uncertainty contour is shaded. Best parameter 
fits are To ~ 110 and g ~ 1.9 for \ 2 ~ 1-9- Both the thick and thin 
shell regimes favour mildly relativistic solutions with £ ~ 1. The 
FS + RS model (lower panel, dashed line) improves the quality 
of the light curve fit by a factor of ~ 2.5 as compared with the 
FS model (lower panel, solid line). Results are consistent with the 
PK02 constraint of To > 68 (arrow) 



Figure 6. GRB 991208: Constraints on To and g assuming a 
thick shell. The best fit value is marked by a solid dot and the 
approximate 1 a uncertainty contour is shaded. Best parameter 
fits are To ~ 130 and g ~ 3.5 for a x 2 ~ 1-5. Both the thick 
and thin shell regimes favour mildly relativistic solutions with 
£ ~ 1. The FS + RS model (lower panel, dashed line) improves 
the quality of the light curve fit by a factor of ~ 2 as compared 
with the FS model (lower panel, solid line). Results are consistent 
with the PK02 constraint of Tq > 68 (arrow). 



tivistic RS with £ ~ 1. At 8.5 GHz, a thin RS is unable 
to provide solutions which improve the FS fit. The best re- 
sult gives Xfs+rs/Xfs ~ 1-4 f° r values of To ~ 110 and 
g ~ 1.5. Thick shell equations fail to converge and therefore 
predict only a lower limit of To > 130. Combining the two 
solution sets should generally enable further constraints. For 
this case, however, combining the sets merely reproduces the 
15 GHz solution since this set of results is statistically dom- 
inant. Therefore, we quote the 15 GHz results as our best 
fits for this RS solution (Figures 5 and 6). Notice that the 
results are consistent with the PK02 limit of To > 68. 



Panaitescu & Kumar 2000) . We find that the fit can be im- 
proved by a factor of ~ 2 by including emission from the 
RS. Both thick and thin shell solutions give Xfs+rs ~ 4.2. 
The thin shell solution favours a mildly relativistic RS with 
best fit values of To ~ 680 and g ~ 1.55 (Figure 7) while the 
thick shell solution gives a highly relativistic solution with 
To ~ 1200 and g ~ 1.5 (Figure 8). We note that this solution 
corresponds to a peak time less than ~ 2 sec which is con- 
siderably smaller than the estimated light width of the shell 
~ 15 s. Both thick and thin regimes produce RS solutions 
which are consistent with the PK02 constraint of To > 150. 



3.5 GRB 991216 

Early radio observations of GRB 991216 suggest the emis- 
sion was already in decline after t ~ 1 day (Frail et al. 2000) . 
As a result, efforts to fit the radio light-curve with a stan- 
dard FS emission model proved difficult. Broadband fits by 
Panaitescu & Kumar (2000) underestimate the early radio 
observations and claim that interstellar scintillation is essen- 
tial in explaining the departures between observations and 
model fluxes. Alternative models have also been suggested, 
including a dual fireball, and a RS flare (Frail et al. 2000, 



3.6 GRB 000418 

Broadband afterglow fits to the FS emission are in agree- 
ment with the radio observations (Berger et al. 2001; PK02). 
It is clear that the addition of a RS contribution will not im- 
prove the fit to the radio light-curve. A minimum flux con- 
tribution is expected from the RS for Po ~ 150 and g ~ 1.5 
assuming a thin shell (Figure 9). This implies a mildly rel- 
ativistic RS with £ ~ 2. The radio flux in the thick shell 
solution decreases with increasing bulk Lorentz factor and 
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Figure 7. GRB 991216: Constraints on To and g assuming a 
thin shell. The best fit value is marked by a solid dot and the 
approximate 1 a uncertainty contour is shaded. RS parameters 
which produce the best fit to the data arc To ~ 680 and g ~ 
1.55 for x 2 ~ 4.2. The lower panel compares the best fit FS+RS 
light curve model (dashed line) and the FS model alone (solid 
line). This value of To is consistent with the PK02 constraint of 
To > 150 (arrow). Note that the thin solution favours a mildly 
rclativistic RS. 



Figure 8. GRB 991216: Constraints on To and g assuming a 
thick shell. The best fit value is marked by a solid dot and the 
approximate 1 a uncertainty contour is shaded. RS parameters 
which produce the best fit to the data are To ~ 1200 and g ~ 1.5 
for a thick shell. The lower panel compares the best fit FS+RS 
light curve model (dashed line) and the FS model alone (solid 
line). This value of To is consistent with the PK02 constraint of 
To > 150 (arrow). In comparison with the thin shell constraint, 
the thick solution tends towards the relativistic limit with £ < 1. 



thus we can only determine a lower limit of To > 50 for this 
regime. It is emphasised that the x 2 ratio never falls below 
1 (with the best fit Xfs+rs ~ 9-7) f° r all values of To and g. 
It also should be noted that the relatively large uncertainty 
region for this burst is due to the predicted faintness of the 
radio flare at the epoch of observations (~ 10 days). 

3.7 Caveats 

A number of caveats apply to our analysis. First and fore- 
most, the hydrodynamical and kinematical constraints pre- 
sented here have been calculated under the assumption that 
the radio flare arises from the reverse shock component of 
the blast wave (MR99 and SP99). Within the grounds of the 
FS model, the radio flare (and in particular those observed 
in 990123, 991216 and 991208) can be accounted for only 
if the characteristic frequencies can be evolved much faster 
than that given by standard dynamics (as argued also by 
Kulkarni et al. 1999 for 990123). In contrast, the RS model 
provides a natural and consistent explanation for the radio 
flare. 

In the framework of "standard" afterglow models - by 
which we mean models that assume a power-law-shocked 



particle spectrum and a constant fraction of energy in elec- 
trons and magnetic fields relative to the thermal energy of 
the shocked particles - we have estimated the RS spectra 
by assuming that the equipartition values for the FS com- 
ponent are unchanged across the contact discontinuity. Al- 
though the validity of this assumption is questionable, it is 
a convenient method of estimating these illusive parameters 
until more is known about the microphysics of GRB shocks. 
We test the strength of this assumption by altering es (the 
parameter for which there is larger uncertainty) by two to 
three orders of magnitude. The effect of strongly increasing 
6b is only a slight shift in the To - g confidence regions. 
Such shifts, however, are always found to lie within the un- 
certainty regions given for each burst, thus supporting the 
validity of our results. 

The premise of our work has been that the physical 
parameters derived by PK02 provide a robust description 
of the FS emission. Most importantly, we assume that the 
PK02 parameters are not significantly biased by neglecting 
the effects of a RS component. This assumption is supported 
by the fact that the PK02 results are weighted towards the 
optical and X-ray data where the sampling frequency was 
large, the observational errors were small, and the RS com- 



© 0000 RAS, MNRAS 000, 000-000 



8 A. M. Soderberg & E. Ramirez- Ruiz 



Table 1. Constraints on To and g obtained by fitting a RS + FS model to early radio observations. The reduced \ 2 f° r the best FS+RS 
model is given as well as the ratio Xps+rs/^fs- For Xps+Rs/-^FS ^ ^ ' s i m P rove d by including emission from the RS. 



GRB 


To 


Thin regime 

g{9low — 9high) XfS + RS 


Xfs+rs/^fs 


To 


Thick regime 

9(9low — 9high) XfS+RS 


Xfs+rs/Xfs 


980519 1 


490 


1.5 


12 


1.4 


>230 


1.5 




- 


990123 


850 


3.4 


0.3 


0.01 


240 


3.3 


0.3 


0.01 


990510 


270 


1.5 


8.5 


8.0 


>230 


1.5 






991208 


110 


1.9 


1.9 


0.5 


130 


3.5 


1.5 


0.4 


991216 


680 


1.5 


4.2 


0.6 


1200 


1.5 


4.2 


0.6 


000418 


150 


1.5 


9.7 


1.0 


>50 


1.5 







ponent was negligible on timescales ~ 1 day. With sparse 
observations and large error bars, the radio data clearly does 
not constrain the multi-frequency fit at the same level as the 
optical. The RS predictions presented here, for all bursts in 
our sample, proved to be generally consistent with the ob- 
served light-curves, and it is critical to note that none of the 
RS predictions significantly over-estimated the observed ra- 
dio emission, even at early times. Our RS predictions there- 
fore provide consistency checks for PK02 physical parame- 
ters, which we find to be in agreement. 



4 THE ELUSIVE F 

As discussed in SR02, the bulk Lorentz factor plays a criti- 
cal role in the majority of GRB emission models (see Piran 
1999 for a review). The Lorentz factor, however, has also 
proven to be one of the most difficult physical parameters 
to constrain from gamma-ray burst observations. Based on 
theoretical predictions and sparse observations, present es- 
timates stretch over two orders of magnitude: r ^10— 10' 5 
(e.g. Meszaros, Laguna & Rees 1993; Wang, Dai & Lu 2000; 
Ramirez- Ruiz & Fenimore 2000; PK02; LS01; SR02). Here 
we present constraints on To for 6 of the best-sampled bursts 
observed to date. Although the Fo constraints are still quite 
broad in some cases (e.g. 000418), information may still be 
gained from the collective Fo distribution for our data set. 
As a method of investigating this, we treat each burst's 
constraint as a probability distribution and produce a to- 
tal histogram by summing over the 6 normalised distribu- 
tions. Specifically, for each burst, we construct a one dimen- 
sional probability histogram for To, given by pi(T)dT. Here, 
Pi(T)dT represents the frequency of each dr value within 
the 1-a To — g confidence region. The pi(r)dr probability 
histograms are each normalised to unity before summing 
them to create a total probability histogram, constructed as 
p(r)dr = J2iPi( T ) dr - Figure 10 plots the total probability 
histogram for the data set as a shaded and smoothed curve. 
A Gaussian fit to the total histogram is over plotted with 
a mean value of To = 150. It should be noted that the dis- 
tribution is weighted heavily to the To constraints derived 
for GRB 991208, since this burst was the most tightly con- 
strained within the data set (due to additional observations 



at 15 GHz). To remove this burst from the total probability 
histogram leaves 5 bursts with radio flare constraints derived 
from 8.46 GHz data only. The total probability histogram 
associated with this smaller set of bursts has a larger mean 
bulk Lorentz factor which may represent the collective con- 
straints more accurately, with To = 270. We emphasise that 
this mean value is consistent with lower limits on Fo as de- 
rived by LS01 for a set of 12 bursts and also with estimates 
for To by PK02 for 10 bursts. As additional To constraints 
are gathered from future bursts, the distribution of Lorentz 
factors should offer vital clues concerning the diversity of 
GRB source expansion velocities. 



5 GRB ENVIRONMENTS AND 
PROGENITORS 

Models for the GRB afterglows indicate that the emission 
comes from a region ~ 10 16 — 10 18 cm from the source of the 
explosion. The nature of the material depends on the GRB 
progenitors, which are at present not known. The question 
of a wind versus a constant density medium is a crucial one, 
since massive stars, one of the leading candidates for GRB 
progenitors (e.g. Woosley 1993; Paczynski 1998), should be 
surrounded by a p oc r~ 2 wind (unless a supernova explosion 
occurs before the burst; Vietri & Stella 1998; Konigl & Gra- 
not 2002; Guetta & Granot 2002). In contrast, GRBs result- 
ing from compact star mergers, the other leading candidate, 
are expected to be surrounded by the interstellar medium 
(e.g. Lattimer & Schramm 1976; Paczynski 1986; Narayan, 
Paczynski & Piran 1992; Perna & Belczynski 2002; Rosswog 

6 Ramirez- Ruiz 2002). 

Considerable discussion has recently been given to 
whether the expected r~ 2 density structure for a stellar wind 
is compatible with analysis of the afterglow light curves (e.g. 
CL00; Lazzati & Perna 2002; PK02; Price et al. 2002). PK02 
found that in half of the cases they modelled, a homogeneous 
ambient medium accommodates the afterglow emission bet- 
ter than the wind-like medium. Although wind-like external 
medium solutions were reported as being a poor match to 
the data in the other cases, PK02 include them as possi- 
ble alternatives for 3 of the 6 bursts in our sample: 991208, 
991216, and 000418. To adopt the r" 2 solutions for these 3 
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Figure 9. GRB 000418: Constraints on To and g assuming a 
thin RS. The best fit value is marked by a solid dot and the 
approximate 1 a uncertainty contour is shaded. A RS contribution 
can be minimised for values of To ~ 150 and g ~ 1.5. Thick shell 
equations provide a lower limit of To > 50 (filled arrow) which 
is less constraining than To > 90 as derived from jet modelling 
(PK02). This shock is mildly relativistic with £ ~ 1.7. The lower 
panel compares the resulting best fit FS+RS model light curve 
(dashed line) with the FS model (solid line). 
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Figure 10. Distribution of GRB Lorentz factors. This distribu- 
tion is essentially a smooth histogram of the data, but one that 
takes into account the uncertainties in the measurements. The 
dashed curve is the distribution under the curve (shaded) but 
smoothed with a Gaussian with a mean value of To = 150. 



bursts would place their RSs in the radiative (fast cooling) 
regime. At best, a wind medium radiative RS can attain an 
optical peak flux which is comparable in magnitude to the 
homogeneous medium adiabatic case. For typical parame- 
ters, however, the optical flux is 6 times fainter in an wind 
environment than the constant density medium case (CLOO). 
Specifically, the peak flux is given by 

where X is the fractional abundance of hydrogen (~ for 
Wolf-Rayet winds) and A is a constant which describes the 
density profile of the wind (CLOO). This peak corresponds 
to an optical flash of magnitude my « 12 as compared with 
an r~ 2 estimated peak of mv ~ 9. Furthermore, the RS 
emission dies off more rapidly when surrounded by a wind 
external medium. The ability to detect a RS contribution 
within early radio observations is then generally reduced 
for an r~ 2 burst. Therefore, after modelling the early af- 
terglows of 991216 and 991208, which are better described 
with the addition of a RS contribution, we argue that the 
homogeneous medium solution is preferred. As a consistency 
check, one can instead assume the lower limits in the value 
of To given by PK02 and fit for the density of the external 
medium. By doing so, we find density values that are con- 
sistent with those from the FS model, which implies that 
there is no direct evidence for large scale clumps within the 
density profile. 

If the progenitors are massive stars then there is an 
analogy to the explosions of core collapse supernovae, for 
which there is abundant evidence that they interact with 
the winds from the progenitor stars. In most supernova 
cases, the radial range that is observed is only out to a 
few 10 17 cm, such that the mass loss characteristics have 
not changed significantly during the time that mass is sup- 
plied to the wind (CLOO). The density in the wind depends 
on the type of progenitor. Red supergiant stars, which are 
thought to be the progenitors of Type II supernovae, have 
slow dense winds. Wolf-Rayet stars, which are believed to 
be the progenitors of Type Ib/c supernovae and possibly 
of GRBs (e.g. MacFadyen & Woosley 1999), have faster 
lower density winds. Deceleration due to this wind starts 
in earnest when about half the initial energy is transferred 
to the shocked matter, i.e. when it has swept up F^" 1 times 
its own rest mass. The typical mass where this happens is 
M oj = E/{T 2 c 2 ) ~ 5 x VOT 6 E 53 Y^%Mq (Meszaros & Rees 
1993). 

Figure 11 shows the (isotropic equivalent) swept up 
mass as a function of radius derived from our r estimates. 
It should be noted, however, that the afterglows sample a re- 
gion < 10 17 cm (~ 0.03 pc) in size. Depending upon the wind 
history of a Wolf-Rayet star during its last few centuries, the 
density structure in this region could be quite complicated 
as the star enters advanced burning stages unlike those in 
any Wolf-Rayet star observed so far (the uncertainty in the 
evolution of massive stars leaves open the possibility of inter- 
action with denser material at early times; see e.g. Ramirez- 
Ruiz et al. 2001). Still, there is no obvious way for the ejected 
mass at 10 17 cm to be much lower than about 1O _5 M0 (as- 
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Table 2. Wind-bubble parameter constraints for 6 bursts. Estimates derived by requiring the expanding observed blast wave to propagate 
in the low-density, (uniform) shocked wind region. Best fit parameters of both M e j and deceleration radius are used (see Figure 11), 
which together with the inferred circumburst density give v w , M and no via equations (14), (15) and (16). 



GRB v w M no 

(1000 km/s) lO- 6 M /yr cm" 3 



980519 
990123 
990510 
991208 
991216 
000418 



0.58 
0.50 
1.4 
1.0 
1.8 
1.4 



0.3 
0.015 
0.2 



100. 
3.4 
55. 



0.7 > 10 4 
0.08 > 10 4 
0.9 > 10 4 



10 



10 



10 



10 




r / pc 



Figure 11. Constraints on the swept up mass M e j as a function 
of radius. The dashed curve is the ejected stellar mass assuming 
M = lO -6 M0yr" 1 and v ~ 10 3 kms _1 . The errors in M e j and 
radius take into account the uncertainties in the Tq estimates. 



suming M = lO" 6 M yr" 1 and v w ~ 10 3 km s _1 ; see dashed 
line in Fig. 11). The low swept up mass inferred is thus 
problematic for the collapsar model. This has lead to the 
suggestion that the fireball expansion may be taking place 
inside the constant density medium that is expected down- 
stream from the termination shock of the massive star wind 
(Wijers 2001; see also Scalo & Wheeler 2001). The radius of 
the wind termination shock at the inner edge of the wind 
bubble can be found by balancing the wind ram pressure 
with the post-shock cavity pressure. The termination shock 
radius, R t , is thus given by 



Rt 



0.4M!r^°n„-j /10 tfpc 



(14) 



where te is the lifetime of the star in Myr and no is the 
interstellar gas density. The density in the uniform shocked 
wind region, n sw , at late times is given by 



3M 



t 



= 0.06M%°n; 



4/5 3/5 6/5.-4/5 



6 "0,3 u w,3 l/ 6 



(15) 



which shows that even if the progenitor star is embedded in 
a dense molecular cloud the observed blast wave can prop- 
agate in a low-density, uniform medium (Wijers 2001). The 
mass within the 1/r 2 wind, Mt, is 



M t ~3x 10" 



,13/10 -9/10 



1 6 



"w,3 



'0,3 



3/w tl /5 M Q . 



(16) 



Comparison with estimates in Figure 11 show that if the 
wind is especially weak (i.e. M ~ lO _6 M0yr _1 ) or the sur- 
rounding pressure is high (no > 10 3 ), Rt falls within the 
range of the relativistic expansion. Considering the radial 
requirements relevant to our sample, we conclude that this 
applies to most cases (see Table 2). Models and observa- 
tions of Galactic Wolf-Rayet stars, however, show that the 
swept-up shell of a red supergiant material at the outer ra- 
dius is at a distance > 3 pc from the star (Garcia-Segura 
ct al. 1996). This radius is sufficiently large that the inter- 
action with the free 1/r wind is expected over the typical 
period of observation of afterglows. The calculations above 
demonstrate that a blast wave expanding into a wind bub- 
ble only works for massive stars with relatively low mass 
winds (< 10 _6 Meyr _1 ) and /or embedded in dense molecu- 
lar clouds no > 10 3 cm~ 3 . 

Although the interstellar and wind models are the two 
main types of environments considered for afterglows, there 
is a different scenario involving a massive star in which the 
supernova explosion occurs before the GRB (e.g. Vietri & 
Stella 1998; Dado et al. 2002; Konigl & Granot 2002; Guetta 
& Granot 2002). The supernova would expand into the pro- 
genitor wind, creating a complex circumburst region in the 
inner part of the wind. Konigl & Granot (2002) have re- 
cently shown, for the case of a pulsar-wind bubble, that the 
shocked wind has a roughly uniform density, similar to that 
found in the normal interstellar medium. 



6 DISCUSSION 

Accepting the inference that the radio flares arise from the 
reverse shock, we place constraints on To, g, and hence on 
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the density profile of the medium in the immediate vicin- 
ity of bursts 980519, 990123, 990510, 991208, 991216, and 
000418. Unlike the continuous forward shock, the hydrody- 
namic evolution of the reverse-shocked ejecta is more sensi- 
tive. As we demonstrated, the temperature of the reverse- 
shocked fluid is found to be non-relativistic for most of these 
bursts. The hydrodynamics of the cold shocked ejecta is very 
different from that of the hot ejecta which is described by 
the BM76 solution. Surprisingly, both cases predict rather 
similar light curves, with decay laws that vary in relatively 
narrow ranges. Hence, it is unlikely that the non-detection 
of radio flares is linked to the relativistic regime of their re- 
verse shock. As argued above, the detection of such flares, 
or firm upper limits, would play an important role in dis- 
criminating between cold and hot shell evolution. 

Moreover, the strong dependence of the peak time of 
this RS emission on the bulk Lorentz factor To provides a 
means to measure this elusive parameter. Indeed, we find 
the fit quality of the FS + RS model to be most sensitive to 
the parameter value of t pea k. Our To estimates are found to 
typically lie between 100 and 10 3 , well above the lower limits 
derived from the requirement that gamma-ray bursts be op- 
tically thin to high-energy photons (see LS01). Constrained 
loosely by equation 1, the peak time is predicted to have a 
lower limit comparable to the light-width of the shell. Un- 
fortunately, the shell width can only be approximated based 
on the duration of the prompt 7-ray emission, and conse- 
quently the lower limit on £ pea k involves a significant degree 
of uncertainty in r . For this sample of bursts, the i pea k val- 
ues are constrained to a large distribution between 1 and 
100 seconds. An increase in the detection efficiency for RS 
flares and flashes will allow t pca k to be better constrained. 
Although large localisation errors plague attempts at early 
(t < 1 day) radio observations, we have shown that signa- 
tures of the radio flare may still be detectable (on top on the 
FS emission) at times significantly after 1 day. Therefore, we 
promote the use of this constraining method in cases where 
a late-time radio flare component would otherwise go unno- 
ticed. Certainly it is true that earlier radio observations are 
the key to improved flare detection, but until current obser- 
vational constraints are lifted, we argue that it is sufficient 
to use observations at t > 1 day and still recover useful RS 
constraints. 

Knowing To, we have then estimated the (isotropic 
equivalent) swept up matter at the radius where the after- 
glow is produced to be < 1O _5 M0 (it should be noted that 
the afterglows sample a region ~ 0.03 pc in size). Provi- 
sional upon the density of the interstellar medium, this may 
be comparable to the wind termination radius; but unless 
the the wind is especially weak or the surrounding pres- 
sure is extraordinary high, the low densities masses inferred 
here may be problematic for the collapsar model. It is evi- 
dent that the environment surrounding a massive star at the 
time of its death is a very opulent one. The complex density 
structure we see in SN 1987A could be a hint of what exists 
in some GRBs progenitors - rich behaviour with multiple 
possible transitions in the observable part of the afterglow. 
If we see such transitions, they can be fairly constraining 
on the properties of the progenitors. So far, the only can- 
didate for having shown a shock transition is the flare-up 
of 970508 one day after the trigger. Due to the deficiency of 



early data, this interpretation is only one of many allowable. 
The prompt optical and radio emission from future bursts 
will give the opportunity of investigating such transitions in 
detail. 

In summary, we show the potential of modelling early 
radio afterglow observations in an attempt to understand 
the physics of the RS. By placing constraints on the prop- 
erties of GRB shocks we hope to enable improved obser- 
vational strategies for future flash and flare detections. We 
show that although current RS theory predicts a detectable 
radio flare component for a wide range of Fo and g param- 
eter values, it appears to be observationally supported for 
only some GRBs. As the sample of well-observed GRB af- 
terglows increases, this method will become applicable to a 
larger set of bursts thereby forming a more statistically sig- 
nificant data set. Analysis of a larger set will offer improved 
diagnostics on both the initial Lorentz factor of the colli- 
mated fireball and also the density of the burst immediate 
environment. 
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